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Abstract. 

We recall that diagonals of rational functions naturally occur in lattice 
statistical mechanics and enumerative combinatorics. In all the examples 
emerging from physics, the minimal linear differential operators annihilating these 
diagonals of rational functions have been shown to actually possess orthogonal 
or symplectic differential Galois groups. In order to understand the emergence 
of such orthogonal or symplectic groups, we analyze exhaustively three sets of 
diagonals of rational functions, corresponding respectively to rational functions 
of three variables, four variables and six variables. We impose the constraints 
that the degree of the denominators in each variable is at most one, and the 
coefficients of the monomials are 0 or dil, so that the analysis can be exhaustive. 
We find the minimal linear differential operators annihilating the diagonals of 
these rational functions of three, four, five and six variables. We find that, even 
for these sets of examples which, at first sight, have no relation with physics, their 
differential Galois groups are always orthogonal or symplectic groups. We discuss 
the conditions on the rational functions such that the operators annihilating their 
diagonals do not correspond to orthogonal or symplectic differential Galois groups, 
but rather to generic special linear groups. 
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1. Introduction 

In previous papers mm it has been shown that the n-fold integrals corresponding 
to the n-particle contributions of the magnetic susceptibility of the Ising model HU 
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[5], as well as various other u-fold integrals of the “Ising class” [Bl|7], or n-fold integrals 
from enumerative combinatorics, like lattice Green functions [iunniiii], correspond 
to a distinguished class of functions generalizing algebraic functions: they are actually 
diagonals of rational functions mmiiiiiiiiis]. As a consequence, the power series 
expansions of the analytic solutions at a; = 0 of these linear differential equations are 
“Derived From Geometry” m and are globally bounded [18) . which means that, after 
a rescaling of the expansion variable, they can be cast into series expansions with 
integer coefficients [HE]. 

In another paper uni revisiting miscellaneous linear differential operators mostly 
associated with lattice Green functions in arbitrary dimensions, but also Galabi-Yau 
operators and order-seven linear differential operators corresponding to exceptional 
differential Galois groups, we showed that these irreducible operators are not 
only globally nilpotent but are also homomorphic to their (formal) adjoints. 

Considering these linear differential operators, or, sometimes, equivalent operatorsU, 
we showed that, either their symmetric square or their exterior square, have a rational 
solution [lOj . This is a general result: an irreducible^ linear differential operator 
homomorphic to its (formal) adjoint is necessarily such that either its symmetric 
sguare, or its exterior square has a rational solution, and this situation corresponds 
to the occurrence of a special differential Galois group [10]. We thus defined the 
notion of being “Special Geometry” for a linear differential operator if it is irreducible, 
globally nilpotent, and such that it is homomorphic to its (formal) adjoint [20]. Since 
many “Derived From Geometry” m n-fold integrals (“Periods”) occurring in physics 
are seen to be diagonals of rational functions mm, we address several examples of 
(minimal order) operators annihilating diagonals of rational functions, and remark 
that they also seem to be, almost systematically, associated with irreducible factors 
homomorphic to their adjoints m- 

Finally in a last paper m revisiting an order-six linear differential operator, 
already introduced in [10] . having a solution which is the diagonal of a rational 
function of three variables, we saw that the corresponding linear differential operator 
is such that its exterior square has a rational solution, indicating that it has a selected 
differential Galois group, and that it is actually homomorphic to its adjoint. We 
obtained the two corresponding intertwiners giving this homomorphism to the adjoint 
relation. We showed that these intertwiners are also homomorphic to their adjoints [21] 
and have a simple decomposition, already underlined in a previous paper [10] , in terms 
of order-two self-adjoint operators. From these results, we deduced a new form of 
decomposition of operators for this selected order-six linear differential operator in 
terms of three order-two self-adjoint operators. We generalized this decomposition 
to decompositions in terms of arbitrary self-adjoint operators of arbitrary orders, 
provided the orders have the same parity [21] . 

This natural emergence in physics of n-fold integrals that are diagonals of rational 
functions, such that their associated linear differential operators correspond to selected 
differential Galois groups, SO{n, C) or Sp{n, C) (or subgroups, like exceptional 
groups [To]), was illustrated on important problems of lattice statistical mechanics 
like the n-fold integrals and of the square Ising model [20], or non-trivial 
lattice Green functions examples m 

The occurrence of diagonals of rational functions necessarily yields linear 

t In the sense of the equivalence of linear differential operators, see CSj. 

i Minimal order operators annihilating diagonals of rational functions may be irreducible, in general. 
In that case, all their irreducible factors have the above property. 
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differential operators that are m globally nilvotenW^. with series-solutions that are 
globally boundec^ [18j . and such that these series identify modulo prime, or even 
powers of prime^ to series expansions of algebraic functions mm- We may call 
these (generally transcendental) holonomic functions “quasi-algebraic” transcendental 
functions. 

If the natural emergence in physics of diagonals of rational functions also 
corresponds to highly selected linear differential operators, the fact, that we observe on 
all our (large number of quite non-trivial) examples of lattice statistical mechanics [H 
[n Hi m m [Ml US], enumerative combinatorics [Mill], that such selected linear 
differential operators have systematically special differential Galois groups, is not 
well understood. The occurrence of selected differential Galois groups {SO{n, C), 
Sp(n, C), or subgroups of SO(n, C), Sp(n, C), like exceptional groups such as G 2 
see m) is now understood, after nni, as the fact that the corresponding linear 
differential operators are non-trivially homomorphic to their adjoints. However the 
“quasi-algebraic” character of the diagonal of rational functions (reduction to algebraic 
functions modulo primes or powers of primes) is not enough to yield the property, for 
the corresponding linear differential operator, to be homomorphic to its adjoint. 

We have accumulated a very large number of example^ of diagonals of 
rational functions for which, quite systematically, the corresponding linear differential 
operators are homomorphic to their adjoints, thus having selected differential Galois 
groups. The set of diagonals of rational functions is an extremely large set: an 
accumulation of examples, that one can hope to be representative of the generic case, 
is just a way to build educated guess, or intuition. Among such results the exhaustive 
analysis of certain sets is worth mentioning. For instance, a set of 210 explicit linear 
differential operators annihilating periods, which are actually diagonals of rational 
functions, arising from mirror symmetric^ (associated with reflexive 4-polytopes 
defining 68 topologically different Calabi-Yau 3-folds, see Ell HU) obtained by 
P. Lairez ES], has been analyzed in El]- Among these 210 operators many correspond 
to “standard” Calabi-Yau DDEs, already analyzed in various papers [7] (I . However, 
remarkably, the other linear differential operators are (non classical Calabi-Yau) higher 
order operators of even orders N = 6, 8, 10, • • • , 24. It was found that all these 
linear differential operators have symplectic differential Galois groups, Sp{N, C), with 
a remarkable canonical decomposition [21) . in terms of self-adjoint order-two operators, 
except the self-adjoint “rightmost operator” which is systematically of order four. 

After these accumulations of examples of diagonals of rational functions 
yielding systematically selected differential Galois groups, one might be led to 
conjecture that the linear differential operators annihilating diagonals of rational 
functions are necessarily homomorphic to their adjoints. Such a conjecture is in 
fact trivially false, as can be seen with the simple 3 F 2 hypergeometric example 

ft And, thus, rational number exponents, for all the singularities of these Fuchsian equations, 
wronskians of these linear differential operators that are A^-th roots of rational functions, etc ... 
i They can be recast into series with integer coefficients mm- 
§ And, in fact, modulo any integer. 

II Unpublished results. 

:( Using a smoothing criterion of Namikawa m, Batyrev and Kreuzer found m 30241 reflexive 4- 
polytopes such that the corresponding Calabi-Yau hypersurfaces are smoothable by a fiat deformation. 
In particular, they found 210 reflexive 4-polytopes defining 68 topologically different Calabi-Yau 3- 
folds with hii = 1. 

f They are order-four irreducible operators satisfying the Calabi-Yau condition: they are, up to a 
conjugation by a function, irreducible order-four self-adjoint operators m- This amounts to saying 
that the exterior square of these order-four operators is of order five. 
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3^2 ([1/3, 1/3, 1/3], [1, 1], 3® a:) which is the Hadamard product [TJ [U [521130] of three 
times the simple algebraic function (1 — 3^ (Hadamard cube), and is, thus, the 

diagonal of a rational function mm, the corresponding order-three linear differential 
operator having a SL(3, C) differential Galois group. In other words this operator, 
associated with the diagonal of a rational function, cannot he homomorphic to its 
adjoint (even with an algebraic extension). One can find many other similar counter¬ 
examples of diagonals of a rational function with a SL{3, C) differential Galois group. 
For instance, 3F2([l/3, 1/3, 1/5], [1, 1], 3'^ 5'^ x) which is the Hadamard product oi two 
times (1 — 3^x)“^/® with (1 — 5^x)“^/®,or 3F2([l/2, 1/2, 1/3], [1, 1], 2'^3'^x) which 
is the Hadamard product of two times the algebraic function (1 — 2^a;)“^/^ with 
(1 — 3^ are two such examples. 

The property, for a linear differential operator, to be homomorphic to its adjoint, 
can be seen to be “some kind’fif] of Poincare duality [31I3I]. Actually, considering 
a diagonal of a rational function, amounts to considering an algebraic variety. If this 
algebraic variety is “smooth enough” (not too singular), then one will have a Poincare 
dualitjlil, which, in some abstract H-module perspective [34], should correspond to 
the previous “homomorphism to the adjoint” property. A simpler occurence of this 
phenomenon is shown by Bogner in proposition 3.4, page 5 of [33j . In simpler words, 
the emergence of a selected differential Galois group should be natural for the diagonals 
of a rational function, when the rational function is not “too singular”. 

For diagonals of rational functions that we know to occur naturally in theoretical 
physicsQ, understanding the emergence of such a duality yielding selected differential 
Galois groups, at least in a physicist’s perspectiv^ requires the analysis of other 
well-defined sets of diagonals of rational functions, along the line of the previously 
mentioned exhaustive analysis mm of the 210 explicit linear differential operators 
arising from mirror symmetries. 

This paper will provide such an analysis, introducing well-defined sets of diagonals 
of rational functions of three, four and six variables, showing that all these examples 
yield selected differential Galois groups, namely orthogonal or symplectic groups. We 
will first analyze a large set of rational functions of three variables, namely rational 
functions with denominators with degree bounded by 1. It turns out that all these 
diagonals of rational functions of three variables correspond to modular forms, so 
that all can, in principle, be written as 2 T 1 hypergeometric functions with two pull¬ 
backs [35]. 

We will finally provide miscellaneous examples that are worth keeping in mind 
when one tries to understand in which case selected differential Galois groups do not 
occur for diagonals of rational functions, so that one could see if such “exceptional” 
cases can also occur in a physics framework. 

ft To be more explicit would require to write a mathematical paper, calling out a Deligne-Steenbrink- 
Zucker theorem saying that Gauss-Manin connections are “variations of polarized mixed Hodge 
structures”, the associated graded modules being Gauss-Manin connections of smooth projective 
varieties, and that they are self-adjoint by Poincare duality. This is far beyond the scope of this 
“learn-by-example” paper (see for instance ED). 

t The Poincare duality theorem is a result on the structure of the homology and cohomology groups 
of manifolds (see chapter 4, page 53 of [32]). 

I And we know why, see mm- 

^ Who wants to understand why orthogonal, symplectic or exceptional groups occur in his problems. 
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Let us recall the definition of the diagonal of a rational function of n variables 
7Z(xi,..., Xn) = V(xi,..., Xn)/Q(xi,..., Xn) where V and Q are polynomials of 
xi, ■ ■ • , Xn with rational coefficients such that Q(0,..., 0) 0. The diagonal of TZ is 

defined from its multi-Taylor expansion 

OO CXD 

n(^Xi, X2, ■ ■ ■ , xjj = X! ■■■ X! , (1) 

mi = 0 m-n — 0 

as the series of one variable 

OO 

Diag 5 ^ 2 ; • ■ - ? Xn^ ^ ^ ^ Rm, ...,m ' X . (2) 

m — 0 

In the generic case, denoting d the degree of the denominator of a rational function 
of n variables, and the order of the minimal order linear differential operator 
annihilating the diagonal of this rational function, the order Q. grows with n and d 
like ~ d” (see e.g. [5S]1: 

n= i. ((d-i)"+i + (- 1 )-+!. (d-i)^ < d". (3) 

From a theory of singularity perspective m. one can imagine that the analysis of 
such diagonals of rational functions depends “essentially” on the denominator of the 
rational function. For that reason we will restrict most of the examples of this paper 
to rational functions with a numerator normalised to 1 , so that the number of cases 
to analyze will be reasonable. Note, however, that, even in a theory of singularity 
perspective, this restriction to numerators equal to a constant, is not innocent, as will 
be seen below (see the Remark in section (ILD). We impose this restriction, for the 
simplicity of the calculations, and in order to perform exhaustive analysis of certain 
sets of examples. 

Furthermore, in order to be able to find a linear differential operator for the 
diagonal, using reasonable computer resources, we will, for a given number of variables 
n, restrict to denominators of the lowest possible degree, imposing, for instance, that 
the degree, in each variable Xi, of the monomials xf^ x^^ ■ ■ ■ x'f" of the denominator 
of the rational function is at most 1. Formula ([3]) gives the upper bound for the order 
of the linear differential operator, which is actually reached when the polynomial 
at the denominator, has all its monomials (no monomial has a zero coefficient). In 
practice, in the denominator of the rational function examples emerging from physics, 
the polynomial at the denominator is (fortunately) sparse, being the sum of a quite 
small set of monomials. In such (physical) cases, the order of the corresponding linear 
differential operator is less than the one given by ([3]) which grows like ~ d" and thus 
becomes quickly too large for any formal calculation. 

Imposing a constraint on the degree in each variable, instead of a constraint 
on the degree of the denominator d = di -I- d 2 -f • • • + d„, reduces the number of 
monomials and, consequently, reduces quite drastically, the order of the corresponding 
linear differential operator. 

Furthermore, it is clear (by definition of the diagonal of a function) that scaling 
the variables Xi ^ Xi ■ Xi, amounts to performing a simple scaling on the diagonal: 
A(a:) —> A(A • x) where A = Ai ■ A 2 • • - An. Consequently we will, in this paper, 
often restrict the coefficients of the monomials to a narrow set of small integer values. 
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namely 0, 1, or 0, ±1, in order to have a reasonably small number of cases to analyze 
to avoid any explosion of the combinatorics. It is, of course, clear that the coefficients 
of rational functions cannot, in general, be reduced to 0, ±1. 

Diagonals of rational functions of two variables necessarily yield algebraic 
functions [n I Hi [sa l40l 1^ . Therefore, in the next sections, we will consider 
diagonal of rational functions, of more than two variables. 


3. Diagonals of rational functions of three variables 


Let us, first, consider diagonals of rational functions of the form 1/(1 — P{x, y, z)) 
where P{x., y, z) is a polynomial of three variables x, y, z, sum of monomials a;™ y" 
where the degrees m, n, p are 0 or 1, and where the coefficients in front of these 
monomials are restricted to take two values 0, 1. We will say that two rational 
functions of the form 1/(1 — P) are in the same class if they have the same diagonal 
and hence the same linear differential operator annihilating this diagonal. With these 
two constraints (on the degrees and values of the coefficients), one finds only 20 
classes of rational functions (or diagonals). For all these cases, the linear differential 
operator is of order two, the diagonal being a 2 F 1 hypergeometric function, or a 
HeunG function |35) that can be rewritten as a 2 F 1 hypergeometric function with two 
possible pullbacks., which means that this is, in fact, a modular form |35) . 

With the same constraint on the degrees m, n, p to be 0 or 1, restricting the 
coefficients in front of these monomials to take three values 0, ±1, one obtains 85 
classes of linear differential operators. 


Let us just give a few examples, the exhaustive list of results being given in a web 
page of supplementary material |42) . 

• For the polynomial P{x, y, z) = x+y+z—xz—yz+xyz, the diagonal of 
1/(1 —P{x, y, z)), which corresponds to the sequence [1, 3, 13, 63, 321, 1683, ■••] of 
Central Delannoy numbers (see Sloane’s on-line encyclopedia [43] of integer sequences: 
A001850), is a simple algebraic function (1 —6a: -I- annihilated by an order- 

one operator. 


• For the polynomial P[x, y, z) = x+y+z—xz—xyz, the diagonal of 
1/(1 — P{x, y, z)) is the pullbacked 2 F 1 hypergeometric function 

(1 - 12a: + 14a:^ -h 12x^ ■ 2 Fi{\^, ^], [1], Pi{x)^ 


where: 

This Hauptmodul 
in |46] 1 

Pi (a:) = 


Pi(x) = 1728- 


(1 — 11 a: — a:^) 


-. (4) 

(1 — 12 a: -I- 14 a:^ + 12 a:^ -|- a:^)^ 

Pi (a:) corresponds to the Hauptmodul 12^/j5 (see Table 5 


123 . ^5 


53 • 


(.»+250*+3125F *■ = 1-llx-i^' 

This diagonal corresponds to the sequence [1, 3,19,147,1251,11253, • • •] of Apery 
numbers (see Sloane’s on-line encyclopedia [33] of integer sequences: A005258). The 
pullbacked 2 P 1 hypergeometric function (jTj) can also be written with another pullbaek 

(1 +2281 +4941= -228i’ + . 2 ^ 1 ( 1 /^, ^1, [1], fliF)) 

where; ftd) = 1728,- id-lli-iT-^ 

(1 -b 228 a:-b 494 a:2 - 228 a:3 -Fa:4)3 
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This Hauptmodul [35] corresponds to the Hauptmodul 12 ^/(see Table 4 

in [33] 1 


P2{x) 


123 . 

(z2 + lOz +5)3 


with: 


53 • X 

1 — 11 a: — 


(7) 


Changing Pi{x) into P 2 {x) amounts to changing z -H- 53/z. 

These two pullbacks, Y = Pi{x) and Z = P 2 {x), are related by an algebraic 
curve, namely the (genus zero) modular curve: 

254 . 53 . . y 6^6 110949 . 2^^ ■ 11^ • 53 ■ ^5^5 . (y + Z) 

+ 3 • 235 • . y4^4 . (2735484611275 + Z'^) - 107937074856652 TZ) 

+ 2^5 ■ 5^ ■ ^3^3 . (Y + z)- (4046657341273198 {Y^ + Z^) 

+ 523793662474799327- YZ) 

+ 3 • 238 ■ 5 ■ y^Z^ ■ (13^ ■ 647451979 ■ {Y"^ + Z^) 

+ 53- 3482348755357972227- Y^Z^ - 16391442082714013450- (yz3 +y3^)) 

+ 30720- yz- (y + Z)- (36983- (Y^ + Z'^) +2421471845930417- (yz3 + y3y) 
+ 506711612589929401008y^Z^) 

+ (y 6 _^2-5) -246683410950- (YZ^ +Y^Z) -441206965512914835246100- y3z3 
+ 383083609779811215375 - (y^Z^ + Z'^y^) 

+ 2^3 - 33 - 5^ - y Z - (y + Z) - (2535689 - {Y^ + Z^) + 134848657695982 yZ) 

- 59719680- yZ- (227547- (y^ + Z^) -83299968230- yZ) 

+ 2^3. 33 °. 5-31-yz-(y + Z) -2^4. 332 -yz = o. (8) 


The pullbacked 2 F 1 hypergeometric function (|T]), or equivalently ([51), is a modular 
form (see the J 5 of Maier [46]). 

• For polynomial P{x, y, z) = x+y+xz + yz+xyz, the diagonal of 
1/(1 — P(x, y, z)) is the pullbacked 2 F 1 hypergeometric function 

(1 - 20 a: +54a:^ -20a:3 + x^)-^/^ ■ 2 F,{[^, ^], [ 1 ], Pi(a:)) 


where: 


Pi (a;) = 1728- 


x^ ■ (1 — 18 a: + a;^) (a; — 1 )^ 


(1 -20a; +54a;2 -20a;3 + a;4)3' 

This Hauptmodul Pi(x) corresponds to the Hauptmodul 123 /j 2 (see Table 5 in (46 

232 . ^2 


(9) 


Pi (a:) = 


123 . ^2 


with: 


2 = 


( 10 ) 


{z + 256)3 “ “ (1 — a;)2 (1 — 18 a; + a;2) 

This diagonal corresponds to the sequence [1,5,49,605,8281,120125---] (see 
Sloane’s on-line encyclopedia of integer sequences: A243945). Similarly, the 
pullbacked 2 F 1 hypergeometric function ([^ can also be written 

(l-20a; + 294a;2-20a;3 + a;4)-3/4. 2Pi([^, ^], [1], P2(x)) 

- (1 — 18 a; + a;3)2 • (1 — x)^ 


where: 


P2{x) = 1728- 


( 11 ) 


(1 -20a; + 294x2 - 20a;3 +x4)3- 
This Hauptmodul P 2 {x) corresponds to the Hauptmodul 123/^2 (see Table 4 in [46]) 
123. z 


P2{x) = 


{z + 16)3 


with: 


(1 — a ;)2 (1 — 18a; + a; 2 )' 


( 12 ) 
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Changing Pi{x) into P 2 {x) amounts to changing z o jz. These two pullbacks, 
Y = Pi(x) and Z = P 2 {x), are related by an algebraic curve, namely the (genus 
zero) modular curve: 

59 y3^3 . 12 y2^2 . {Y Y Z) - Y Z ■ (16 - 4027 -YZ +16 Z'^) 

-6A- {Y + Z)- {Y^+ IA&1YZ + Z"^) +2^^ ■ 2,^ ■ YZ = 0. (13) 

The pullbacked 2 F 1 hypergeometric function ([S]), or equivalently is a modular 
form (see the J 2 of Maier |46], or even Ja of Maier |46] but with z = • (1 — 

x)V(l -18x +x2)2). 


• For the polynomial P(x, y, z) = xyz + xy + xz + yz + x + y + z, the diagonal 
of 1/(1 — P(x, y, z)) is the pullbackec0 2 P 1 hypergeometric function 


1 


2 Pi([^, ^], [1], P(a:)) where: P(x) = 


1 — X 


54 X 


\3 • 


(14) 


(l-x)- 

This diagonal corresponds to the sequence [1,13,409,16081,699121, 32193253 ••• ] 
(see Sloane’s on-line encyclopedia of integer sequences: A126086). This series can 
also be written as the pullbacked 2 P 1 hypergeometric function 

(1 (1 a - a 29 x + 9 x ^ 2i^i([Y2> [il- a(2^)) 

a(x) = 3456. »’^(1 -Vi+3,t“-,t4)» 


where: 2 . 

(1 — x)^ (1 -|-429x -|-3x^ — x^)^ 
or as the pullbacked 2 Pi hypergeometric function 

(1 -x)-^/^- (1 -51x +3x^ -x^)-^/"^- 2 Pi([^, ^], [1], P 2 {x 


(15) 


where: 


P 2 (x) = 13824- 


(1 — 57x -I- 3x^ — x^) 


(16) 


(1 — x)9 (1 — 51X -I- 3 x^ — x^)^' 

These two pullbacks correspond respectively to the Hauptmodul [3S] 12^ jjs of Table 
4 in [?5] and 12 ^/j 3 of Table 5 in [45] : 

123 • z 193 . ^3 

= 7—— 7 —n7v7> -^ 2 ( 3 ;) = 


(z -k27) • (z -k3)3’ 

with: z = 2 ■ 


(z -k 27) • (z -k 243)3 


93 


(17) 


1 — 57x -I- 3x^ — x3 

Changing Pi(x) into P 2 (x) amounts to changing z 3^/z. Again, these two 
pullbacks are related by a (genus zero) modular curve: 

227 . 5 ® • y32’3 . (y + z) + 2^3 • 5® ■ • (27 Y^ - 45946 PZ + 27 Z^) 

-f 2® • 53 • 3® • YZ - {Y + Z)- (y3 + 241433 PZ -k Z^) 

-f 729 ■ {Y'^ + Z*) - 779997924 • (yz3) -h Y^Z + 31949606 ■ ■ Y'^Z'^ 

+ 2^ ■ 3^^ ■ 31- Y Z ■ {Y + Z) - 2^2 . 3^^ ■YZ= 0. 

The pullbacked 2 Pi hypergeometric function (BH) is a modular form (see the J 3 of 
Maier |46]b 

• For the polynomial P(x, y, z) = xyz + xy + xz + yz, the diagonal of 
1/(1 — P(x, y, z)) is the pullbacked 2 P 1 hypergeometric function 


1 


2 Pi([^, ^], [1], P(ai)) where: P(x) = 


1 — X 


27 x2 
(l-x)3- 


(18) 


% Note that this pullback can be obtained using the Maple program ”hypergeomdegS” of V. J. 
Kunwar and M. van Hoeij |44l I45| . 
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This diagonal corresponds to the sequence [1,1, 7, 25,151, 751 • • ■ ] (see Sloane’s on-line 
encyclopedia of integer sequences: A208425). 

In [Appendix A| we give a set of diagonals of 1/(1 — P{x, y, z)) which can be 
seen to be modular forms: they can be written as 2 F 1 hypergeometric functions with 
two different pullbacks. 


Most of the other examples of diagonals, given in a web page as well as in 
the supplementary material |42| . are not yet in Sloane’s on-line encyclopedia of 
integer sequences. For instance the sequence [1, 4, 42, 520, 7090, 102144, • • •], which 
corresponds to the diagonal of the rational function 1/(1 — x — y — xy — x z — y z), 
can be written as a 2 F 1 hypergeometric function with a rational pullback: 


1 


. 1 — 16 a; — 


/r 

) - 2^1 ([ 


1 5 

12 ’ 12 


I, [1], 1728- 


• (l-bx)2(2 - 34x- 27x2)' 
(1 — 16 X — 8x^)3 ^ 


All these results show a close relation between diagonals of our simple examples 
of rational functions of three variables, and modular forms. This suggests, quite 
naturally, to perform similar calculations, but, now, with four variables. 


4. Diagonals of rational functions of four variables 

Let us now consider diagonal of rational functions of the form 1/(1 — P(x, y, z, w)) 
where P(x, y, z, w) is a polynomial of four variables x, y, z, w, sum of monomials 
y'^ z^ where the degrees m, n, p, q are 0 or 1, and where the coefficients in 
front of these monomials are restricted to take only two values 0, 1. With these 
two constraints (on the degrees and on the values of the coefficients), one finds an 
exhaustive list of only 879 casefl 

Remark 1: We have used the package “HolonomicFunctions” written by 
C. Koutschan 071 SHI SS], based on the method of creative telescoping 001 EH, 
which enables to obtain directly, and very efficiently, the linear differential operator 
annihilating the diagonal of a given rational function, without calculating the series 
expansion of the corresponding diagonaQ. From time to time Koutschan’s algorithms 
do not provide the minimal order linear differential operator. It is thus necessary 
to systematically check whether the linear differential operator obtained by creative 
telescoping is minimal, and if not, to find the minimal one. 

Remark 2: The command “FindCreativeTelescoping”, described in EH, is 
(usually) extremely fast, but it uses some heuristics, which means that sometimes 
it can return a non-minimal result or run forever: it is not an “algorithm” in the strict 
sensei- 


Among these 879 cases, the linear differential operators annihilating the diagonals 
have minimal orders running from I to 10, as given in Table [T] 

% When different rational functions 1/(1 — y, z, w)) yield the same diagonal (in practice the 
same first ten coefficients of the series) we select one rational function to represent the diagonal. This 
is the way we define these 879 different classes of rational functions. 

I One does not obtain the linear differential operator from a “guessing procedure” on the series of 
the diagonal. Furthermore this algorithm is “certified”; we are sure that the operator annihilates the 
diagonal. 

f In contrast, Chyzak’s algorithm |52| is designed such that it finds the minimal order operator, 
but it is often much slower than Koutschan’s heuristics (in particular for the very large operators 
emerging in physics H [3 ESI [M]). In Koutschan’s package, Chyzak’s algorithm is implemented in 
the command CreativeTelescoping (same input/output specification as FindCreativeTelescoping). 
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Table 1. Number of operators corresponding to the various orders. 


Order of Oper. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Number of Oper. 

1 

2 

20 

128 

240 

231 

155 

54 

41 

7 


The order-one linear differential operator corresponds to the diagonal of the 
rational function 1/(1 — z — wy — xy). It is easily seen that its diagonal is equal 
to 1 , so the order-one operator is just Dx- 

For the order-two linear differential operators, the corresponding diagonals are 
two pullbacked 2^1 hypergeometric functions, respectively 

(1 -32x +16 x2)-i/4- 2 Fi([^, A], [1], P{x)) 

-3. (2 -71a; +Wx^) 


where: 


P{x) = 1728■ 


(1 - 32 a; -f 16a;2)3 ’ 

for the diagonal of 1/(1 — {y + z + x z + xw + xyw)) which reads 

1 -f 8 a; -f 156a;2 -f 3800a;3 -f 102340 a;'‘ -f 2919168 a;® -f 86427264 a;® 
-f 2626557648 a;^-f 81380484900 X®-f 2559296511200 a;® -f •• 


and 


(1 -40x -f 16x®) ^], [1], P(x)) 

• (1 — 44x — 16 X®) 


where: 


P{x) = 6912- 


(1 — 40x -I- 16x2)3 

foiQ the diagonal of 1/(1 — (y -f z -f xw + xzw + xyw)) which reads: 


1 -flOx -f 246x® -f7540x® -f 255430 x^ 
-b 13042646760 x^ -b 508236930630 x® - 


-b 9163980 X® -b 341237820 X® 
■ 20101587623260X® -b ••• 


(19) 

( 20 ) 
( 21 ) 


( 22 ) 


Note, however, that the series (ISHI) of the diagonal of 1/(1 — (y-bz-bxz-bxrc-b 
xyw) ) actually identifies with the diagonal of a rational function of just three variables 
1/(1 — X — y — z — xy), already found among the previous 20 cases of section ([3]). 
Similarly the series (l22l) of the diagonal of 1/(1 — (y -bz + xw + x zw -bxyw)) also 
identifies with the diagonal of a rational function of three variables 1/(1 — x — y — 
z - xy -yz). 

The results for all these 879 cases are given exhaustively in our web page of 
supplementary material |42) . Let us summarize these results in the following. 

For all the twenty cases, corresponding to order-three linear differential operators, 
we have SO{3, C) differential Galois groups. As a consequence [TU], all these linear 
differential operators are actually symmetric squares of order-two operators, some 
with very simple 2 ^^! hypergeometric functions, namely 2 Ai([ 3 / 8 ,1/8], [1], 256x®) or 
2 Fi([ 1 / 3 , 1/6], [1], 108X®), some with, at first sight, more involved HeunG function 
solutions [35] which turn out to be pullbacked 2 F 1 hypergeometric functions, with 
two possible pullbacks, and, in fact, modular forms |35j . 

The 128 order-four linear differential operators are (non-trivially) homomorphic 
to their adjoints. They have 50(4, C) differential Galois groups and have a canonical 

t These two series iiM and are not in Sloane’s on-line encyclopedia http://oeis.org 
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decomposition m of the form {Ai i ?3 + 1) ■ r(a;), where Ai and are respectively 
order-one and order-three self adjoint operators, r{x) being a rational function. 

Similarly, the 240 order-five linear differential operators are (non-trivially) 
homomorphic to their adjoints. They have SO{b, C) differential Galois groups and 
have a canonical decomposition [21] of the form {Ai Bi C 3 -\-Ai Cs)- r{x), where Ai, 
Bi and C 3 are respectively two order-one and one order-three self adjoint operators, 
r{x) being a rational function. 

The 231 order-six linear differential operators are (non-trivially) homomorphic 
to their adjoints. They have SO( 6 , C) differential Galois groups and have a canonical 
decomposition [21] of the form {Ai Bi Ci D 3 -\- Ai Bi Ai D 3 -\- Ci D 3 1) ■ r{x), 
where the Ai, Bi, Ci, operators are order-one self-adjoint operators, the rightmost 
operator D 3 being an order-three self-adjoint operator. 

The following 155 order-seven linear differential operators are (non-trivially) 
homomorphic to their adjoints. They have SO(7, C) differential Galois groups 
and have a canonical decomposition [5T] of the form {Ai Bi Ci Di E 3 -\- Ai Bi E 3 -\- 
Ai D\ E 3 Ai Bi Cl -\-Ci Di E 3 -\- Ai -\- Ci -\- E 3 ) ■ r(^x), where the Ai , Bi , Ci , Di 
operators are order-one self-adjoint operators, the rightmost operator E 3 being an 
order-three self-adjoint operator. 

The 54 order-eight linear differential operators are (non-trivially) homomorphic 
to their adjoints. They have SO{ 8 , C) differential Galois groups and have a canonical 
decomposition described in |21] . generalization of the previous ones, with, again, five 
order-one self-adjoint operators, and a rightmost order-three self-adjoint operator. 

The 41 order-nine linear differential operators are (non-trivially) homomorphic 
to their adjoints. They have SO{9, C) differential Galois groups and have a canonical 
decomposition described in |21] . generalization of the previous ones, with, six order- 
one self-adjoint operators, and a rightmost order-three self-adjoint operator. 

Finally, the seven order-ten linear differential operators are (non-trivially) 
homomorphic to their adjoints. They have iS'O(10, C) differential Galois groups and 
have a canonical decomposition described in [21j . generalization of the previous ones, 
with, seven order-one self-adjoint operators, and a rightmost order-three self-adjoint 
operator. 

These results are reminiscent of the results obtained on a set of 210 explicit 
linear differential operators annihilating diagonals of rational functions, arising 
from mirror symmetries and corresponding to reflexive 4-polytopes [28] . recalled in 
the introduction. One notes, however, that the symplectic Spin, C) differential 
Galois groups with a canonical decomposition in order-two self-adjoint operators 
and a rightmost order-four self-adjoint operator, encountered with these reflexive 4- 
polytopes examples, is now replaced by orthogonal SOin, C) differential Galois groups 
with a canonical decomposition in order-one self-adjoint operators and a rightmost 
order-three self-adjoint operator. 

The calculations performed here, in order to see that these (quite large) linear 
differential operators are (non trivially) homomorphic to their adjoints and to find their 
canonical decompositions [21] , are similar to the ones described in [21] for the reflexive 
4-polytopes examples: in order to find the intertwiner, we introduced a specialized 
algorithm because the Maple command Homomorphisms{adjointiL), L) never 
terminates on these large operators [2T]. We use a fuchsian linear differential system 
associated to L, the theta-system |21j . a slight generalization of the companion system 
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which has simple poles at each finite singularity. One then finds a rational solution of 
an associated system with similar coefficients (its second symmetric/exterior power), 
which gives the intertwiner. An inversion of this intertwiner modulo L gives the 
intertwiner corresponding to the Maple command Homomorphisms{L, adjoint(L)). 
Finally, one obtains the canonical decomposition from simple euclidean divisions |21) . 

For the last seven order-ten operators even the theta-system calculations were 
quite massive, and required up to two weeks of CPU time and up to 80 Gigaoctets of 
memory for one linear differential operator. 

5. Diagonals of rational functions associated with orthogonal as well as 
symplectic groups 

It is tempting to simply, and straightforwardly, generaliz^J these two sets of results 
for diagonals of well-defined finite sets of rational functions, namely the symplectic 
Sp{n, C) differential Galois groups [55] (with an order-four rightmost self-adjoint 
operator as for reflexive 4-polytopes [55]), and the orthogonal SO{n, C) differential 
Galois groups (with a rightmost order-three self-adjoint operator). The situation can 
be slightly more involved (and richer) than a straightforward generalization of the 
previous results. In fact, diagonals are not necessarily solutions of an irreducible 
linear differential operator. 

To see this let us consider the diagonal of the rational function of four variables 
with a factorized denominator R = 1/(1 — x — y — u — z)/{l — u — z — uz), which 
reads: 

Diag(R) = 1 -k 42 a; -k 4878 + 748020 -k 130916310 -k 24762428460 x® 

-f 4929691760532 X® -f 1017691904736936 x’’ -f ■ • • (23) 

It is solution of an order-seven linear differential operator which factorize^ into an 
order-three and an order-four operator: L-j = ■ L 4 . Note that L 4 does not 

annihilate this diagonal. 

This means that the series ([23|) is solution of T 7 = L 3 • L 4 , but not of L 4 which 
has a solution analytic at x = 0 , with a series expansion with integer coefficient^ 
different from 1)23^ : 

1 -f 214x -k 97278x2 + 53983020x3 + 32898451110 x'‘ + 21172639875156x3 
+ 14121624413802444x3 + ••• (24) 

The series L 4 {Diag{R)), solution of L 3 , is globally bounded [HE]- If one normalizes 
L 4 to be an operator with polynomial coefficients (instead of being a monic operator), 
the series L 4 {Diag(R)), solution of L 3 , is a series with integer coefficients: 

17888 + 25769200 X + 17312032256x2 +8722773606816x3 + 3775743401539200 x'^ 
+ 1486619414765913792x3+548416028673746513280x3 + ••• 

f In particular after a set of other unpublished results we have obtained. 

% Note that this diagonal is factorized but is not the Hadamard product of two diagonals (see 
[Appendix B.2| l. 

§ Using the DFactorLCLM command one sees that one does not have a direct sum factorization, just 
the simple factorization = L 3 • L 4 . 

ji A natural question corresponds to ask if such series dSU with integer coefficients are necessarily a 
diagonal of rational function. It is true for ODEs of minimum weight for the monodromy filtration m, 
but in the general case, it is still a conjecture. It, however, seems that one can prove that such series 
are “automatic” (i.e. reduce to algebraic functions) modulo powers of primes (G. Christol, private 
communication). 
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Note that the order-three operator L 3 has a canonical SO{3, C) orthogonal 
decomposition [U] {AiBiCi + Ai + Ci) ■ r{x), where Ai, Bi, Ci, are order- 
one self-adjoint operators. In contrast L 4 has the canonical S'p(4, C) symplectic 
decomposition {A 2 B 2 + 1 ) • p{x), where A 2 , B 2 are order-two self-adjoint operators. 
In other words, the diagonal (|23)) is associated with both orthogonal and symplectic 
differential Galois groups! This is, in fact, the situation we expect generically: the 
diagonal of the rational function will be solution of a non-irreducible linear differential 
operator, each of its factors corresponding to orthogonal or symplectic differential 
Galois groups. 

We have, for instance, the same results for another rational function such that its 
denominator is factorised. The diagonal of R = 1/{1 —x — y — z — u)/(l — u — z —xz), 
namely 

1 +44x +5061x^ + 771000 x^ -h 134309890 -h 25316919264 

-H 5026804760628 X® -h • • • (25) 

is solution of an order-seven linear differential operator which factorizes into an order- 
three and an order-four operator: T 7 = T 3 • L 4 (but, again there is no direct sum 
factorization). This diagonal (1^51) is solution of L^ but not of L 4 . If one normalizes L 4 
to be an operator with polynomial coefficients (instead of being a unitary operator), 
the series L 4 {Diag{R)), solution of L 3 , is a series with integer coefficients: 

16 - 94464 X - 127052100x2 -86146838400x3 -44244836836200x4 

- 19495756524980736x3 - 7791904441995369696 X® -h • • • (26) 

The order-three linear differential operator L 3 is MUM (maximal unipotent 
monodromy [55]) and has a canonical SO{3, C) decomposition {Ai Bi Ci -\- Ai -\- 
Ci) ■ r(x), where Ai, Bi, Ci, are order-one self-adjoint operators. It is not the 
symmetric square of an order-two operator, it is homomorphic to the symmetric square 
of an order-two operator. In contrast L 4 has the canonical 5*^(4, C) decomposition 
{A 2 B 2 + 1 ) • p{x) where A 2 , B 2 are order-two self-adjoint operators. 


6. Diagonals of rational functions corresponding to nFn-i hypergeometric 
functions 


Beyond these finite sets of examples of diagonals of rational functions, namely the 
210 reflexive 4-polytopes operators [35] with symplectic Sp{n, C) differential Galois 
groups, and these 879 operators of section ([3]) associated with diagonals of rational 
functions of four variables with orthogonal SO{n, C) differential Galois groups, one 
can find infinite families of diagonals of rational functions for which exact results can 
be obtained corresponding to nFn-i hypergeometric functions. 

For instance, the diagonal of the rational function of three variables R = 
1/(1 — X — 0 — y"), with n being a positive integer, is a 2 n^2ra-i hypergeometric 
function 


^*“3(7- - -= 

\l—X — z — 

2 n 2 n ^ 2 n + V 2 n + l’ 2 n + l" ’2n+V’ 

1 1 2 2 3 3 n-1 n-1 n (2n-hl)2"+4 

L 5 5 5 5 7 7 ''‘7 7 7 J7 O ^ 

n n n n n n n n n 
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The corresponding linear differential operator is an order-2n linear differential operator 
having a symplectic Sp{2 n, C) differential Galois group. 

The diagonal of the rational functions of three variables R = 1/(1—x—z —x" y”), 
with n being a positive integer, is a nFn-i hypergeometric function 



(28) 


The corresponding linear differential operator is an order-n linear differential operator 
having an orthogonal 0(n, C) differential Galois group. 

The diagonal of the rational functions of three variables R = 1/(1 — x — z—xy"), 
with n being a positive integer, is a hypergeometric function. For n even it 

reads: 



(29) 


The corresponding linear differential operator is an order-n linear differential operator 
having a symplectic Sp(n, C) differential Galois group. 

For n odd the same formula (0^)) holds. Note, however, that the argument 1/2 
appears in both the first and second list of arguments of the nFn-i hypergeometric 
function, and, hence, can be avoided, which is thus a n-iFn -2 hypergeometric 
function. The corresponding linear differential operator is an order-(n — 1) operator 
having a symplectic Spiri — 1, C) differential Galois group. 

7. Diagonals of rational functions associated -with operators non 
homomorphic to their adjoints 

After all this accumulation of examples of diagonals of rational functions associated 
with orthogonal, or symplectic, differential Galois groups, it is tempting to 
conjecture that counter-example^ like 3F2([l/3, 1/3, 1/3], [1, 1], 3®x), are quite 
“rare exceptional cases” that one can easily detect, and, hopefully, understand, 
as situations of algebraic varieties “sufficiently singular” to break the Poincare 
duality [32]. Along this line of seeking for diagonals of rational functions with a 
“sufficiently singular” denominator, let us try to provide examples of operators that 
are not homomorphic to their adjoints. 

Appendix B provides two attempts to find examples of diagonals of rational 
functions such that their corresponding linear differential operators would not 
be homomorphic to their adjoints. We first study in [Appendix B.l| a set of 
singular denominators for the rational functions, the polynomials of three variables 
corresponding to classifications of singular varieties performed by V. L Arnold |56j . 

§ Given in the introduction corresponding to a SL(S, C) differential Galois group, and thus the 
corresponding operator cannot be homomorphic to its adjoint even with an algebraic extension. 
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All the corresponding operators yield symplectic Sp{n, C) differential Galois groups: 
the kind of singular behaviour required to “break the Poincare duality”, considered in 
V. 1. Arnold [S0] for three variables, is not “sufficiently singular”. 

It could be that these “sufficiently singular” situations require much more than 
three variable^ Therefore we have considered, here, rational functions of sia; variables, 
and since the theory of singularities of algebraic varieties [SzlET] suggests that the 
situation of algebraic varieties “singular enough” often correspond to denominators 
that factoJttl. we have considered, in [Appendix B.2[ rational functions of the form 
1/((1 — Pi)(l — P 2 )) where Pi and P 2 are two polynomials of two different sets 
of three variable^. Again we found that aZ|| the operators that could be analyzed 
among the 170 operators in [Appendix B.^ are homomorphic to their adjoints, and 
have symplectic or orthogonal differential Galois groups. 

We also considered the (minimal order) linear differential operators annihilating 
large sett0 of Hadamard products [H [H [29] of algebraic functions^ and, each time., 
obtained that these linear differential operators are homomorphic to their adjoints. 
Surprisingly, finding diagonals of rational functions such that their annihilating 
linear differential operators are not homomorphic to their adjoints, is not so easy, 
examples like 3P2([l/3, 1/2, 1/2], [1, 1], 12^x), being quite rare. Let us revisit these 
simple hypergeometric examples yielding operators that are not homomorphic to their 
adjoints. 


7.1. Examples of diagonals of rational functions with no homomorphism to the 
adjoint 


The hypergeometric function 3P2([l/3, 1/2, 1/2], [1, 1], 12^a;), actually corresponds 
to a SL{3, C) differential Galois group. It can be seen as the diagonal of a rational 
function of six variables 


jj. ( _1 -9xy _\ 

V(l — By — 2a; + 3?/^ + 9a;2j/)- (1 — u — z) ■ (1 —v — w)) 

= 1 + 12 X + 648 + 50400 + 4630500 x"^ + 468087984 x*^ 


(30) 


The corresponding (order-three) linear differential operator is not homomorphic to 
its adjoint, even with an algebraic extension. This diagonal (I30|) is of the form 
Diag{Ri{x, y) ■ R 2 {u, z) ■ Rsiv, w)), (where i?i, R 2 and R 3 are simple rational 
functions). It is, thus, the Hadamard product [T] H] [23] of the three diagonals 
Diag{Ri{x, y)), Diag{R 2 {u, z)) and Diag{R 2 {v, w)), which are simple algebraic 
functions, respectively: 

*®^(l _ 2 X — 3y + 3y2 + 9x2 y) (1 — 9x)i/3 ’ ^ 


t Counter-examples, like 3F2([l/3, 1/3, 1/3], [1, 1], 3®x), correspond to diagonal of algebraic 
functions of three variables (which are very simple since they are products of the algebraic functions 
of one variable), but this means that they are diagonals of rational functions of, at first sight, six 
variables. 

tfThis is actually the case for the counter-examples, like 3i^2([l/3, 1/3, 1/3], [1, 1], 3® ai) mentioned 
in the introduction. 

I This also corresponds to considering the Hadamard product of the diagonal of 1/(1 — -Pi(a:, y, z)) 
and of the diagonal of 1/(1 — P 2 {u, v, w)). 

i Even using Koutschan’s creative telescoping program m, two operators among the 170 were 
difficult to obtain. Among these 170 operators, 71 are so large that they cannot be analyzed using 
the DEtools commands, and required to switch to a theta-system approach (see m)- 
^ Some miscellaneous heuristic examples are given in the supplementary materials m- 









Diagonals of rational functions 


16 


Diag(- - ') = Diag(- -) = --—— 7 ^. (32) 

^Vl —u — zJ ^Vl —V —wJ (1 — 

The linear differential operators, annihilating the diagonals (EB and (l32l) . are 
very simple order-one linear differential operators. There is, of course, no relation 
between their differential Galois groups and the SL{3, C) differential Galois group 
of the order-three operator annihilating their Hadamard product (I30p : the Hadamard 
product does not preserve algebraic structures, like the differential Galois group. 


Remark: Let us consider the following rational function of six variables: 
_ N{x, y, z, u, V, w) _ 

{1 — 2x — 3y -1-9-1-3y^) - (1 —z — u) ■ (1 — v — w) 

1 — 9 xy 1 1 


(33) 


-b 


1 — 2x — 3y -1-9x^y -1-3y^ ' 1 — z — u ' 1 — v — w 

It is straightforward to see that the diagonal of this rational function (1551) of six 
variables is nothing but the sum of the three diagonals ED) and El, namely: 

1 2 
(1 -9x)i/3 + (1 -4x)i/2- 


This diagonal is solution of an order-two linear differential operator which is 
(obviously) a direct sum: 





(35) 


The differential Galois group of the (order-three) operator that annihilates (1501) . 
and the one of the (order-two) operator (1551) that annihilates (1551) . are quite 
different, even if the denominators of El and El are the same. In a theory 
of singularity perspective 133 El], the statement that the differential Galois group 
depends essentially on the denominatoiQ of the rational function that “encodes the 
singularities”, has to be taken “cum grano salis”. The restriction we have imposed 
on our rational functions by imposing the numerators to be equal to one, is far from 
innocent. 


1.2. More examples of diagonals of rational functions with no homomorphism to the 
adjoint 


A slight modification of the previous example (150|) amounts to considering the diagonal 
of a rational function of six variables 

^((1 — 3y — 2x-|-3y^-|-8x^y)-(I—M — z)-(I—n — w))’ ^ ^ 

where the coefficient of x^ y has been modified (9 x^ y changed into 8 x^ y). The 
corresponding (order-five) linear differential operator is not homomorphic to its 
adjoint, even with an algebraic extension, and its differential Galois group is SL{5, C). 

Another slight modification of (1501) amounts to considering the diagonal of a 
rational function of five variables 


Diag(^ 


1 — 9 X y 


(1 — 3y — 2x-|-3y^-|-9x^y)-(I—M—x — z — w) 


(37) 


I Mathematicians would say that changing the numerator may change the n-form one integrates, 
that one must consider the Gauss-Manin Picard-Vessiot module of these differentials. We just try 
here, heuristically, to make crystal clear that numerators matter. 
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The corresponding (order-four) linear differential operator is not homomorphic to its 
adjoint, even with an algebraic extension, and its differential Galois group is S'T(4, C). 

More examples of diagonals of rational functions with linear differential operators, 
that are not homomorphic to their adjoints, are given in [Appendix C| 

To sum-up: The cases such that the diagonals of rational functions do not 
yield linear differential operators homomorphic to their adjoints, are far from being 
understood, either from a differential algebra viewpoint, or from a theory of singularity 
perspective. From an experimental mathematics perspective, what we see is that 
diagonals of rational functions “almost systematically” (but not always !) yield linear 
differential operators homomorphic to their adjoints, thus giving selected differential 
Galois groups. Diagonals of rational functions were seen mm to naturally emerge 
in physics. However, in a mathematical framework, not related to physics, diagonals 
of rational functions seem to “almost systematically” yield orthogonal or symplectic 
groups, with, at first sight, no obvious “physical interpretation”. 

8. Conclusion 

We have introduced well-defined sets of diagonals of rational functions of three, four 
and six variables, showing that all these examples yield selected differential Galois 
groups, namely orthogonal and symplectic groups. 

It has been seen that, in our set of rational functions of three variables with 
denominators with degree bounded by 1, all diagonals correspond to modular forms 
that can all be written as 2 F 1 hypergeometric functions with two pull-backs |35) . 
related by a modular curve. 

We have seen that a set of 879 diagonal of rational functions of four variables 
correspond to orthogonal SO{n, C) differential Galois groups with a remarkable 
canonical decomposition m with a rightmost self-adjoint operator of order three. 
These results were obtained using the very powerful package “HolonomicFunctions” 
written by C. Koutschan gzniiiin], based on the method of creative telescoping |50j . 
which enables to obtain directly, and very efficiently, the linear differential operator 
annihilating a given diagonal of a rational function, without calculating the series 
expansion of the corresponding diagonal. In order to find the homomorphisms of 
these operators to their adjoints, which is the first step towards the analysis of the 
differential Galois groups of these operators and their “canonical decompositions” [21] , 
we have, for large linear differential operators, also used a new algorithm that requires 
to work on the linear theta-system associated with the operators m- 

We have also seen ( [Appendix B.2| below) that 170 diagonals of rational 
functions of six variables of the form 1/ ((1 — Pi{x, y, z)){i — P^iu, v, w))), actually 
correspond to a quite rich set of linear differential operators with orthogonal or 
symplectic differential Galois groups. The systematic analysis performed in this paper 
of these three sets of diagonals of rational functions of respectively three, four and six 
variables, suggests that diagonals of rational functions “almost systematically” yield 
orthogonal or symplectic differential Galois groups. 

A contrario, we have provided in section © a few miscellaneous examples of 
diagonals of rational functions where the corresponding linear differential operators 
are not homomorphic to their adjoints. These cases, such that the diagonals of rational 
functions do not yield linear differential operators homomorphic to their adjoints, are 
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far from being fully understood. Is it possible that such cases could also emerge with 
the diagonals of rational functions appearing in physics? This remains an open and 
challenging question. 
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Appendix A. Diagonals of rational functions of three variables: some 
modular forms 


Let us recall the two Hauptmoduls [3 
(z +27)- (z +3)3’ 


123 /j 3 and 123/^4 of Table 4 of Maier [46] : 
123- z-(z + 16) 


H4iz) = 


(z^ + 16z + 16)3 


• The diagonal of 1/(1 — x — y — y z — x z) as well as the diagonal of 
1/(1 — X — y — z + xz), correspond to the sequence [1,4,36,400,4900,63504, • • •] 
of the complete elliptic integral A(4a;3/2) (oeis number A002894 in Sloane’s on-line 
encyclopedia) can be written as a pullbacked 2 F 1 hypergeometric function: 

(1 - 16a: + 16^], [1], 774 ( 2 :)) with: z = - —(A.l) 

Recalling Maier’s paper [4^ one knows that (jA.l|) can alternatively also be written 
a pullbacked 27^1 hypergeometric function with the 123/j4 of Table 5 of Maier |46j . 
The diagonal of 1/(1 — x — y — z x z) is clearly a modular form. As a byproduct 
this suggests that the diagonal of the rational function of three variables 

7- ^ , , - — -- (A.2) 

(1 — X — y — z + a;z)' (1 — x — y — xz — yz) 

is zero. This can be checked directly. 

• The diagonal of 1/(1 — x — y — z — xyz), which corresponds to the 
sequence [1, 7,115, 2371, 54091,1307377, • • • ] (oeis number A081798 in Sloane’s on-line 
encyclopedia) can be written as a pullbacked 27^1 hypergeometric function: 

(1 _27a: +30:^ - x^)-^/^ ■ {I - • 27^1 ([^, ^], [1], 773 ( 2 )) 


with: 


or more simply as 
1 




1 — a; 


1 — 30 a: + 3 a;^ — a:3 


21X 


(1 — a:)3 


(A.3) 


(A.4) 
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• The diagonal of 1/(1 — x — y — z — xy — x z — y z), which corresponds to the 
sequence [1,12, 366,13800, 574650, 25335072, • • • ] (no oeis number) can be written as 
a pullbacked 2 F 1 hypergeometric function: 


(1 -48a; 


24a;2)-i/4. 

with: 



(A.5) 


• The diagonal of 1/(1 —x—y —x z—y z +x y z), which corresponds to the sequence 
[1,3, 25, 243, 2601, 29403, • • • ] (oeis number A245925 in Sloane’s on-line encyclopedia) 
can be written as a pullbacked 2 F 1 hypergeometric function: 

(1 - 12a; - lOa;^ - 12a;^ + ^], [1], Hi{z)^ 

with: z = (A.6) 


• The diagonal of 1/(1 — x — y — z + xyz), which corresponds to the 
sequenc^ [1, 5, 67,1109, 20251,391355, • • ■ ] (oeis number A124435 in Sloane’s on-line 
encyclopedia) can be written as a pullbacked 27^1 hypergeometric function: 


(1-ha;) ■ {1 -21x+Zx^ +x^) • 2 -F 1 ([^, ^], [1], 773 ( 2 )) 


with: 


1 — 24 a; -h 3 -h 


(A.7) 


• The diagonal of 1/(1 —x-|-j/-|-z-|-x?/-|-xz — j/z-hxj/z), which corresponds 
to the sequence 

[1,11,325,11711,465601,19590491, ••• ] (no oeis number) can be written as a 
pullbacked 27^1 hypergeometric function: 


(1 _46x +x^)-^/^ • (1 +x)-i/^ • 2 T’i([^, ^], [1], H,{z)) 


with: 


1 — 52x -h x^ 

2 X 


(A.8) 


Appendix B. Two attempts to break the Poincare duality 

Along the line, sketched in section ([7]), which amounts to seeking for diagonals 
of rational functions with a “sufficiently singular” denominator, let us try to find 
examples of linear differential operators that are not homomorphic to their adjoints. 
We study here two sets of singular denominators for the rational functions, first 
polynomials corresponding to classifications of singular varieties performed by V. I. 
Arnold [56], then denominators that factor into two polynomials. 

Appendix B.l. Diagonals of rational functions associated with singular algebraic 
varieties 

If one believes that the situations where the linear differential operator annihilating 
the diagonal of a rational function should correspond to situations of algebraic 

X Number of effective multiple alignments of three equal-length sequences. 
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varieties “singular enough” to break the Poincare duality [31] , it is tempting to study 
the singular algebraic varieties classified by V. I. Arnold [SSj. The classification 
of the simplest singularities turned out to be related to Lie, Coxeter and Weyl 
groups, An, Dn, En, and to the classification of platonic solids in Euclidean three 
spaces (56). Arnold’s paper gives a set of polynomials of three variables that are 
obvious candidates to be considered as denominators of rational functions of three 
variables: P = x'^ z +y^ + +2y z''‘, z + y z'^ + ay^ z, + y^ + z'^ + ax y z'^, 

With these polynomials, we obtained the linear differential operators annihilating 
the diagonals of the rational functions of the form 1/P. These linear differential 
operators of various orders (order 22 for Qio in [56], order 14 for Q 12 in [56], order 
12 for ^ii, order 8 for S' 12 , order 10 for U 12 , order 16 for Z 13 , order 18 for IT 12 , •■■) 
are all homomorphic to their adjoints, their differential Galois groups being Sp{n, C) 
symplectic groups, the canonical decomposition m being in terms of only order-two 
self-adjoint operators. 

Appendix B.2. Diagonals of rational functions of six variables 

Apparently the kind of singular behaviour required to “break the Poincare duality”, 
considered in V. I. Arnold |56j for three variables, is not “sufficiently singular”. It could 
be that these “sufficiently singular” situations require much more than three variables. 
Therefore we have considered in this section, rational functions of six variables 
(see section cu and ( [Appendix C[ ) below). Since the theory of singularities of 
algebraic varieties suggests that the situation of “singular enough” algebraic varieties 
often correspond to algebraic varieties that factor, we have studied exhaustively the 
diagonals of rational functions of the form 1/ ((1 — Pi{x, y, z))(l — P 2 {u, v, w))), 
where the degree of the two polynomials in each of their three variables is less than 
one, and their coefficients are 0 or 1. Note that this also corresponds to considering 
the Hadamard product of the diagonal of 1/(1 — Pi{x, y, z)) and of the diagonal of 
1/(1 — P 2 {u, V, w)). The number of such classes yielding different diagonals is only 
170. We have obtained all the corresponding linear differential operators using the 
Mathematica “HolonomicFunctions” package [47]. The order of the (minimal order) 
linear differential operators runs from 2 to 12. The number of linear differential 
operators corresponding to the various orders is given in Table m- 


Table Bl. Number of operators corresponding to the various orders, number of 
operators with symplectic and orthogonal differential Galois groups. 


Order of Oper. 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Number of Oper. 

2 

3 

19 

13 

39 

0 

52 

0 

36 

0 

6 

Number of Sp(n, C) 

2 

0 

7 

0 

35 

0 

52 

0 

36 

0 

6 

Number of SO{n,<C) 

0 

3 

12 

13 

4 

0 

0 

0 

0 

0 

0 


We also give the number of linear differential operators having respectively a 
symplectic and an orthogonal differential Galois group. Among these 170 linear 
differential operators, 91 can be analyzed using the DEtools command in order to 
see that they are homomorphic to their adjoint, find their differential Galois group 
and their canonical decomposition |21] . The other linear differential operators are 
too large to be analyzed that way: they reguire to switch to a differential theta- 
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system m in order to find the intertwiners to their adjoint, and, then the canonical 
decomposition from a simple euclidean division | 21 j . Note that these calculations 
{theta-system calculations) are still quite massive. 

We found that the three operators with a SO{3, C) orthogonal differential Galois 
group were self-adjoint. 

We found that all the seven order-4 linear differential operators with a symplectic 
differential Galois group (namely 5*^(4, C)) are self-adjoint, or conjugated to their 
adjoints (by a simple rational function), except one with a (L 2 • M 2 -\- 1 ) • r{x) 
canonical decomposition m- We found that the 12 linear differential operators with 
a SO{4:, C) orthogonal differential Galois group had a (Li ■ M 3 -|- 1) ■ r{x) canonical 
decomposition Ell- 

All the 13 order-5 linear differential operators which have an SO{5, C) orthogonal 
differential Galois group have a (Li • Mi ■ N 3 -\- Li -\- N 3 ) ■ r{x) canonical 
decomposition |21| . The four order -6 linear differential operators with an 50(6, C) 
orthogonal differential Galois group have a (Li • Mi • A^i • P 3 -|- Li • P 3 -|- A^i • P 3 -|- 
Pi • Ml -I- 1) • r{x) canonical decomposition |21) . 

We found that the other 35 order -6 linear differential operators have a 5p(6, C) 
symplectic differential Galois group. Among these 35 symplectic operators only two 
have a (P 2 • M 2 ■ N 2 -\- L 2 -\- N 4 ) ■ r{x) canonical decomposition [ 21 ], all the other 
having a (P 2 • M 4 -|- 1 ) ■ r{x) canonical decomposition pT| . 

We found that all the 52 order -8 linear differential operators have a 5p(8, C) 
symplectic differential Galois group with a (P 2 ■ M 2 ■ -\- L 2 -\- N 4 ) ■ r{x) canonical 

decomposition EH- 

All the 36 order-10 and 6 order-12 linear differential operators have respectively 
5p(10, C) and Sp{12, C) symplectic differential Galois groups, with (L 2 • M 2 ■ N 2 ■ 
P 4 + ■ • •) • r{x) and (P 2 ■ M 2 • A ^2 • P 2 ■ Q 4 + • ■ •) • 'f'{x) canonical decomposition [21] . 

Remark 1: In all the decomposition we have obtained, the rightmost self-adjoint 
operator EH, is always of order three (for orthogonal differential Galois groups) or of 
order four (for symplectic differential Galois groups), except an order-4 operatoir with 
a (P 2 • M 2 -I- 1) • r{x) decomposition and the two order -6 linear differential operators 
corresponding to the diagonals of the two rational functions 

{1 — xy — xz — yz) • (1 — v — w — uv — uw), (B.l) 

(1 — z — xy — xz — yz — xyz) ■ (1 — u — v — w — uvw), (B.2) 

the corresponding operators having a (L 2 ■ M 2 ■ N 2 -\- L 2 -\- N 2 ) ■ r{x) canonical 
decomposition EH. 

Remark 2: Using Koutschan’s creative telescoping program [47], three linear 
differential operators were quite difficult to obtain compared to the others. They 
correspond to the following denominators of the rational functions: 

(1 — X — y — z — xz — xyz) ■ (1 — u — v — w — uw — vw — uvw), (B.3) 

(1 — X — y — z — xz — xyz) ■ (1 — u — v — w — uv — uw — vw — uvw), (B.4) 

(1 — X — y — z — xz — xyz) ■ (1 — u — v — w — uv — uw — vw). (B.5) 

Koutschan’s creative telescoping program gives an order-12 linear differential operator 
for dHU). Note that one always needs to verify that the operator obtained from this 
program is actually the minimal order operator annihilating the diagonal^. Actually, 

t It may not be irreducible: see, for instance, the L 7 = L 3 • L 4 operator associated with the diagonal 
(|23]l in section 
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performing Hadamard products muni modulo primes, one can forecast that minimal 
order for (IB.3|) . We found an order 10 which suggests that the order-12 linear 
differential operator obtained by Koutschan’s creative telescoping program |47j , is not 
the minimal order linear differential operatoiQ: this order-12 actually factorises. It is 
the product of two order-one operators and of the minimal order-10 linear differential 
operator annihilating the diagonal: L 12 = Li ■ Mi ■ Liq. The factorization of 
such very large linear differential operatortQ is a quite difficult task. In a first step 
calculations modulo prime have been performed that enable us to find the order of 
the minimal order operator annihilating all these diagonals. We have, however, been 
able to perform these factorizations (in characteristic zero, not modulo primes). 

Remark 3: When the linear differential operator annihilating the diagonal, 
obtained from Koutschan’s program m. is not the minimal order linear differential 
operator, one can try to obtain that minimal order operator by factorization (DFactor 
in Maple). Unfortunately the factorization of (very) large operators like many of these 
operators, cannot be obtained using straightforwardly DFactor DEtools command in 
Maple. The factorization of the largest operators has been obtained, in the following 
way: one first obtains a large set of coefhcients of the series of the corresponding 
diagonals (that is the non-trivial part of the calculation), and then use it as the input 
of a “guessing procedur^^’ . 

Remark 4: Koutschan’s creative telescoping program m gave us an order-9 
linear differential operator corresponding to the diagonal of 1/((1 — x — y — z — x z — 
xy z)(l — u — V — w — uw — uvw)) (i.e. the Hadamard square of the diagonal of 
1/(1 — X — y — z — X z — xy z)). In fact this order-9 linear differential operator is not 
minimal, the minimal order operator being of order 8. This program gave us also four 
order-11 linear differential operators, like, for instance, the operator corresponding to 
the diagonal of (IB.41) . Similarly this order-11 linear differential operator is not minimal, 
the minimal order operator being of order 10 (the degree of the polynomial coefficients 
is 74). This is also the case for the three other, at first sight, order-11 linear differential 
operators corresponding to the following denominators of the rational functions: 

{1 — z — xz — yz — xy) ■ (1 — u — v — w — uw — vw — uvw), (B.6) 

(1 — y — z — xz — xy — yz) ■ (1 — u — v — w — uv — uw — vw — uvw), (B.7) 
(1 — X — y — z — xz — xyz) ■ (1 — u — v — w — uw — vw). (B.8) 

They are not minimal order operators, the minimal order linear differential 
operators, annihilating the corresponding diagonals being of order 10 (the degree of 
the polynomial coefficients being respectively 42, 51, 51). Therefore one finds that 
there is no order-9 or order-11 operators for this set of 170 diagonals. 

The linear differential operator annihilating the diagonal of (IR5I) was the most 
difficult to obtain using Koutschan’s creative telescoping prograirijj. This program 
gives an order-13 linear differential operator. From this exact differential operator one 
can obtain as long as necessary series of the diagonal, and study, modulo some primes, 

I One can even forecast the degree of the polynomial coefficients of that order-10 operator: the 
degree is 74. 

J Using the DFactor command of DEtools in Maple, or any other method. 

% Essentially of the same type as in gfun. 

J The program computes the “telescoper” (the ODE) and the “certificate”. It is a general observation 
that in most examples, the certificate is much larger than the telescoper m- We just need the 
telescoper. An algorithm that computes telescopers without computing the corresponding certificates, 
has already been built [581159) . 
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the linear differential operators annihilating these series. One finds, that way, that the 
minimal order operator annihilating the diagonal of (IB.51) is of order 12 (with degree 
85). We actually factorised this order-13 operator, and obtained the minimal order-12 
operator that annihilates this diagonaQ. 

Remark 5: Seeking for the (large ...) linear differential operators annihilating 
the diagonals of (IB.41) . (IB.51) . it seems natural to take into account the selected form 
of all these diagonals of rational functions of six variables, which are, actually, the 
Hadamard product [HIT] of two diagonals of rational functions of three variables. In 
the case of (IB.41) . (IB.5|) (and even (IB.dll '), the linear differential operators annihilating 
these various diagonals of rational functions of three variables, are simple order-two 
operators. One can imagine to obtain these linear differential operators using the gfun 
command “hadamard product” of two simple order-two operator^. Unfortunately, 
the linear differential operators obtained that way are of order much larger than the 
one obtained from Koutschan’s creative telescoping of a diagonal of a function of six 
variables ! The same remarks apply for the analysis of all the 170 diagonals of Table 

(EH). 


Let us give, here, miscellaneous, simple examples of these diagonal of rational 
functions of six variables (all the exhaustive results being given in a supplementary 
material [H]). 


Appendix B.2.1. Order-two operators: pullbacked 2 F 1 hypergeometric function 
• The diagonal of the rational function of six variables 



— w — uv — uvw) ■ (1 — z — xy)J 

1 +6x +78x^ -b 1260x3 -b 22470x4 -b 424116x3 -b 8305836x® -b 


(B.9) 


is annihilated by an order-two linear differential operator. This diagonal is the 
Hadamard product EUS] of the diagonal of 1/(1 — z — xy) and of the diagonal 
of 1/(1 — ru — uv — uvw), namely the two simple algebraic functions: 

(1-4x)-4/2 (B.IO) 

= 1 -b2x -b 6x2 -b 20x3 -b 70x4 -b 252x3 -b 924x3 . 

(1 -6x -bx2)-4/2 (B.ll) 

= 1 -b3x -bl3x2 4-63x3 _|_ 321 a;4-b 1683x3 -b 8989x3 _|_ ... 


This diagonal is actually a pullbacked 2 F 1 hypergeometric function: 
(1 -24x + 48 x 2 )- 4 / 4 . 2 Fi([^, ^], [l],Pl(x)) 


where: 


Pi(x) 


~^4x -bl6x2) 
(1 -24 X -b 48x2)3 


(B.12) 


% After completion of this work we were told by C. Koutschan that giving extra options to 
“FindCreativeTelescoping” enables the program to find directly the minimal order-12 operator in 
84 hours CPU time, instead of the 116 hours CPU time we used to obtain the order-13 operator. 
With these extra options the telescoper is 232560 bytes when the certificate requires 32743760 bytes, 
to be compared with 265432 bytes and 38422496 bytes for respectively the telescoper and certificate 
in the order-13 calculation. 

I The gfun[hadamardproduct] (eql, eq2, y(z)) command determines the linear differential equation 
satisfied by the Hadamard product of two holonomic functions, solutions of the linear differential 
equations eql and eq2. 
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This Hauptmodul can be seen to be of the form 12^ /(see Table 4 of Maier |46) 1: 


Pi{x) = 12^ 


with: 


1 — 24 a; + 16x^ 


(B.13) 


{z + 16)3 2x^ 

This diagonal (iBJl) has another pullbacked ^Fi hypergeometric function representa¬ 
tion, corresponding to the other pullback P 2 {x) = 12^ • {{z -I- 256)3; 

1 5 


(1 -24a;+528a;2)-i/4- 2Fi{[z^. ^2^, [1], P2{x)) 

- (1 — 24 a; -I- 16 a;^)^ 


where: 


P2ix) = 3456- 


(B.14) 


(1 -24 a; -k 528a;2)3 ’ 

The pullbacked 2 P 1 hypergeometric function (IB.12I1 (or equivalently (IB.14I1 1 is thus a 
modular form. 

Note that this diagonal (IB.9I) is, in fact, actually identical to the diagonal of the 
rational function of three variables 1/(1 — x — y — z — xy + x z). 

• The diagonal of the rational function of six variables 

1 


Diag(^ 


(B.15) 


(1 — w — uv — uvw) ■ (1 — z — xy — xy z) 

= 1 -k 9 X -k 169 -k 3969 x^ -k 103041 x^ -k 2832489 x^ + 80802121 x® -k • • • 

is annihilated by an order-two linear differential operator. This diagonal is the 
Hadamard square of the diagonal of 1/(1 — z — xy — xyz), namely the Hadamard 
square of the algebraic function (|B.llll . This diagonal is actually a pullbacked 2 F 1 
hypergeometric function; 

(1 - 36x + 134 3:^ - 363:^ + [1], Pi ( 3 :)^ 

Pi(x) = 27648- 


where: 


(1 — 36 X -I-134 x2 — 36 x3 -|- x"‘)3 
It is also be written as pullbacked 2 F 1 hypergeometric function: 

(1 +444x + 134x2 +444x3 + 23 Fi([^, ^], [1], P^ 


(B.16) 


where: P 2 {x) = 3456- 


X • (1 — x)2 • (1 — 34x + x2) 


244 


(B.17) 


(1 + 444 X + 134 x2 + 444 x3 + x'*)3 

which shows that this diagonal is a modular form: these two pullbacks can be seen as 
the Hauptmoduls 123 /_) 4 , jf'^ in Table 4 and 5 of Maier [46] : 

z • (z + 16) , x2 — 34 X + 1 


Pi(x) = 123 


where: 


z = 


(B.18) 


(z2+16z+16)3 ■ 2x 

Note that this diagonal (IB.15|) is, in fact, actually identical to the diagonal of the 
rational function of three variables 1/(1 — x — y — z — xy + x z — y z — xyz). 


Appendix B.2.2. Order-four and order-eight operators 
• The diagonal of the rational function of six variables 

Diag(- -= 1 +36x2 +6300 x^ + ... 

V(1 — w — uv) ■ (1 — xy — xz — y z)J 

is annihilated by a linear differential operator of order four, L 4 , which has 
an orthogonal differential Galois group 5'0(4, C), with a simple canonical 
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decomposition EH, L, = ■ L 3 + 12 , where L 3 is an order-three self-adjoint 

linear differential operator: 

L3 = -1) ■ Dl +Sx ■ {86Ax^ -1) ■ Dl 

-h (2868x2 - 1) • + 276x (B.20) 

• The diagonal of the rational function 

1 \ 


Diag(^- 


^(l — uv — uw — vw) ■ (1 — xy — X z — y z)J 

= 1 -h 36x2 -h 8100 -h 2822400 X® -h • ■ • (B.21) 

is annihilated by a linear differential operator of order four, L 4 , which is self-adjoint, 
its differential Galois group being the symplectic group 5'p(4,C). 

• The diagonal of the rational function 


Diagi - 


1 


V(1 — uv — uw — vw — uvw) ■ (1 — xy — x z — y z)> 

= 1 -b 42x2 -b 13590 x"^ -b 7410480 X® -b • • • (B.22) 

is annihilated by a linear differential operator of order eight, Lg, which 
has a symplectic differential Galois group Sp{8^ C), with a simple canonical 
decomposition EH, ^8 = {L2 ■ M2 ■ N4 + L2 + N4) ■ r(x), where L2 and M2 

are two order-two self-adjoint operators and N 4 is an order-four self-adjoint linear 
differential operator. 


Appendix C. Miscellaneous examples of diagonals of rational functions 
with operators that are not homomorphic to their adjoints 

Another slight modification of ((3(111 amounts to considering the diagonal of a rational 
function of six variables 

1 — 9xy 


DiagL 


(C 11 

.(1 — 3y — 2x-b3?/2-b9x^?/)- (1 —z — u -uz)- (1 — v —w)J 
The corresponding (order-four) linear differential operator is not homomorphic to its 
adjoint, even with an algebraic extension, and its differential Galois group is S'T(4, C). 

Similarly, the diagonal of the five variable rational function 

^^((1 — X -\-8y — 27 xy^ —27xy'^ —9xy -b3y^) - (1 — u — v — uz — v z)j 
is annihilated by an order-three linear differential operator which is not homomorphic 
to its adjoint (even with an algebraic extension), and its differential Galois group is 
SL{ 8 , C). 

Similarly, the diagonal of the five variable rational function 

1 — 9 X y 


Diag(^- 


()■ 


.(1 — 3y — 2x-b3y2-b9x^y)-(l — u — u — rc)^’ 
is annihilated by an order-three linear differential operator which is not homomorphic 
to its adjoint, and its differential Galois group is SL{ 8 , C). 

Similarly, the diagonal of the five variable rational function 

1 — 9 X y 


Diagi - 


(C.3) 


V(l — 3y — 2x-b3y^-b8x^y)- (1 — u — v — w)J’ 
is annihilated by an order-five linear differential operator which is not homomorphic 
to its adjoint, and its differential Galois group is SL{5, C). 
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